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ABSTRACT 

Context. Electric dipole emission arising from rapidly rotating Polycyclic Aromatic Hydrocarbons (PAHs) is often invoked to explain 
the anomalous microwave emission. This assignation is based on i) an observed tight correlation between the mid-IR emission of PAHs 
and the anomalous microwave emission; and ii) a good agreement between models of spinning dust and the broadband anomalous 
microwave emission spectrum. So far often detected at large scale in the diffuse interstellar medium, the anomalous microwave 
emission has recently been studied in detail in well-known dense molecular clouds with the help of Planck data. 
Aims. While much attention has been given to the physics of spinning dust emission, the impact of varying local physical conditions 
has not yet been considered in detail. Our aim is to study the emerging spinning dust emission from interstellar clouds with realistic 
physical conditions and radiative transfer. 

Methods. We use the DustEM code to describe the extinction and IR emission of all dust populations. The spinning dust emission is 
obtained with SpDust, which we have coupled to DustEM. We carry out full radiative transfer simulations and carefully estimate the 
local gas state as a function of position within interstellar clouds. 

Results. We show that the spinning dust emission is sensitive to the abundances of the major ions (Hn, Cn) and we propose a 
simple scheme to estimate these abundances. We also investigate the effect of changing the cosmic-ray rate. In dense media, where 
radiative transfer is mandatory to estimate the temperature of the grains, we show that the relationship between the spinning and 
mid-IR emissivities of PAHs is no longer linear and that the spinning dust emission may actually be strong at the centre of clouds 
where the mid-IR PAH emission is weak. These results provide new ways to trace grain growth from diffuse to dense medium and 
will be useful for the analysis of anomalous microwave emission at the scale of interstellar clouds. 

Key words. ISM: general - ISMxlouds - radiative transfer - dust, extinction - evolution 



1. Introduction 

Discovered in the nineties, the anomalous microwave emission 
(AME) has aroused great interest (Kogut et al. 1996 ; Leitch et al. 
|1997 1. First because it appears in a frequency window that is op- 
timal for the detection of the Cosmic Microwave Background 
(CMB) fluctuations. |Draine & Lazariiii] (|T998 ) proposed that 
AME could be caused by electric dipole emission of rapidly ro- 
tating grains: the spinning dust emission. This mechanism is now 
most often invoked to explain the AME and several models have 
been published ( Ali-Hafmoud et al7||2009t |Ysard & Verstraete 



|20lO}|Hoang et al.|2010||Silsbee et al.|201 1 



. The study of spin- 
ning dust could help in understanding the life cycle of interstellar 
dust grains because it may be a new tracer of the smallest grains, 
the interstellar Polycyclic Aromatic Hydrocarbons (PAHs). 

The preference of spinning dust models over other mech- 
anisms is based on several arguments. First, AME is corre- 
lated with dust IR emission and this correlation is particularly 
tight for the mid-IR emission of small grains. Second, AME is 
weakly polarized as expected for PAHs because these grains are 
not supposed to be aligned with the interstellar magnetic field 
2006 Casassus et al. 2008 Lopez-Caraballo 



(Battistelli et al. 



et al. 201 1 1. Third, the shape and the intensity of the AME can be 
reproduced with spinning dust spectra (e.g. Watson et al. 2005; 
Planck Collaboration et al. 201 lb, and many other references). 
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However, the spinning dust emission depends on the local physi- 
cal conditions (gas ionisation state and radiation field) and on the 
size distribution of small grains. Recent observations of interstel- 
lar clouds point out dissimilar morphologies in the mid-IR and in 
the microwave range that may be explained by local variations of 
the environmental conditions (Casassus et al.|200 6 2008; Ysard 



[eTaLl[2010l |Castellanos~et~aT. 



2011 Vidal et al. 20111. In this 



work we study the spinning dust emission of interstellar clouds 
including a treatment of the gas state and radiative transfer. In 
this context we reexamine the relationship between the AME 
and the dust IR emission. 

The paper is organised as follows. In Section [2] we describe 
the models. In Section [3] we detail our method to estimate the 
gas properties (ionisation state and temperature). In Section [4] 
we present the variations of the spinning dust spectrum with the 
gas density and with the intensity of the radiation field. We also 
consider variations of the cosmic -ray ionisation rate as suggested 
by recent observations. In Section[5]we present the spinning dust 
emission with radiative transfer modelling. Finally, we present in 
Section|6]our conlusions. 



2. Models 



Current models of spi nning dust (|Draine & Lazarian|1998||Ali 



Haim oud et~aT||2009| |Ysard & Verstraete||2010| |Hoang etaL 



2010 Silsbee et al. 2011 1 take into account a number of pro- 



cesses for the rotational excitation and damping of the grains: the 
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emission of IR photons, the collisions with neutral and ionised 
gas particles (H i, H n, and C n), the plasma drag (H n and C n), 
the photoelectric effect, and the formation of H2 molecules at 
the surface of the grains. The publically available SpDusQcode 
( [Ali-Haimoud et al.|20"09"l|SiIsbee et al.|201 l) includes the most 
recent developments regarding the gas-grain interactions and the 
grain dynamics (rotation around non-principal axis of inertia). 
The results of SpDust agree well with other models that include a 
more detailed treatment of the IR emission or of the gas-grain in- 
teractions (Ysard & Verstraete 2010 Hoang et al. 2Q1Q| >. SpDust 
is fast and well-suited for coupling to other codes, especially ra- 
diative transfer codes. In the following, we use SpDust to model 
the spinning dust emission. 

In order to estimate dust emission from the mid-IR to the 
microwave range in a consistent way, we coupled SpDust with 
the dust emission model described by Compiegne et al. ( 201 IK 
DustElV0 DustEM is based on the formalism of IPesert et al l 
( |1990| l and includes three dust types: interstellar PAHs, amour- 
phous carbonaceous grains, and amourphous silicates. We used 
the dust populations defined by Compiegne et al. (2011 1 for the 
diffuse, high galactic latitude interstellar medium (DHGL). For 
PAHs we assumed a log-normal size distribution with centroid 
ao = 0.64 nm and width cr = 0.4, with a dust-to-gas mass ratio 
Mp AH IM H = 7.8 x 10~ 4 . 

In current models, the smallest grains (PAHs) carry the spin- 
ning dust emission that is sensitive to the gas density and the ra- 
diation field intensity Go, but also to the ionisation state (abun- 
dance of the H 11 and C 11 ions, noted xh and xq respectively). 
Radiative transfer calculations are performed with the CRlj^Jtool 
Puvela & Pado an 2003; Juvela 2005), to which we have cou- 
pled DustEM and SpDust. CRT is only used to estimate the dust 
temperature and the resulting dust emission from the mid-IR to 
the microwave range. Our treatment of the gas properties is pre- 
sented in Section[3] 



3. Gas state 

As discussed above, the dynamics of spinning dust grains in- 
volves gas-grain interactions and radiative processes. The spin- 
ning dust emission is therefore sensitive to the gas density («#) 
and temperature (T gas ), and to the intensity of the UV radiation 
field traced by the factor Go- In particular the gas-grain interac- 
tions depend on the gas ionisation state, i.e., the abundance of the 
major charged species (electrons, H n, C n, ect.), which primarily 
depends on r gas , and Go but also on the chemistry occurring 
locally. Realistic modelling consequently requires a consistent 
treatment of the spinning motion of the grains and gas ionisation 
state. For the present work, where we consider the influence of 
radiative transfer on the spinning dust emission (see Section |5J, 
we treat the gas ionisation state with a simplified scheme that we 
present below. Using this scheme, we then discuss the influence 
of riH and Go, and look at the effect of an enhanced cosmic -ray 
ionisation rate as suggested by recent observations (see Section 
141 



Available at http://www.tapir.caltech.edu/yacine/ 



spdust/spdust .html 



Available at http : //www . ias . u-psud . fr/DUSTEM 

3 Go is a scaling factor for the radiation field integrated between 6 
and 13.6 eV. The standard radiatio n field corresponds to Go = 1 and to 
an intensity of 1 . 6 x 1 3 erg/s/cm 2 (|Mathis et al.|l983b. 

4 Available at http://wiki .helsinki . fi/display/~mjuvela@ 
|helsinkifi/CRT| 



When the radiation field intensity is low (Go < 1), inelas- 
tic collisions with neutral and ionised species of the interstel- 
lar gas become the dominant processes for the excitation and 
the damping of the grain rotation ( | Ali-Ha imoud et al. 2009; 
|Ysard & Verstraete|2010| >. The ion fractions xh — bhu/")? <md 
Xc = ncn/riH where = «(H 1) + n(H 11) + 2«(H2) accordingly 
need to be carefully determined to perform a quantitative study 
of the variations of spinning dust emission with environmen- 
tal properties. Where CO has not formed (unshielded regions in 
which most of the gas phase carbon is in the form of C 11 or C 1), 
we estimate the electron and ion fractions (x e - n e /n#, xh, and 
xc) by simultaneously solving the H 1 / H n and C 1 / C 11 equilib- 
ria, including the recombination of carbon with H2 (Rollig et al. 
|2006[ ). Furthermore, we take into account the recombination of 
C n with PAHs as described in |Wolfire etaLlpOOBl . 

In neutral gas and neglecting the contribution of helium, the 
ionisation balance of hydrogen including H2 reads 



(Hi,H 2 ) + CR <± (Hn,H+) + e 
Ccr{1-Xh) = x H x e n H a H , 



(1) 
(2) 



where (cr is the cosmic-ray ionisation rate per second and per 
proton and a H = 3.5 X lQ- 12 (r/300K )-°- 75 cm 3 /s is the Hn re- 
combination rate (Rollig et al. 2006). Unless otherwise stated, 



we assume (cr — 5 X 10 _i/ s _1 H _i . In regions where CO has 
not formed, we assume that C 11 is the dominant ionised heavy 
element and write the electron fraction as x e xh + Xc- The C 11 
abundance thus becomes 



1 



1 + x e n H a H l( C R 



(3) 



On the other hand, xc can be derived from the ionisation balance 
of carbon where we take into account the following reactions: 



Ci + hv 
C 11 + e~ 
C 11 + PAH" /PAH 
C 11 + H 2 



Cn + e" (4) 

Ci (5) 

C 1 + PAH/PAH+ (6) 

CH 2 + + hv or CH + + H 1. (7) 



The rate coefficients k r and k x are taken from Wolfire et al 
(|2008|l. Fro m the database of the Meudon PDR code ( |Le Petit 
|et al.|2006| ), we take 



1.7 10 15 + 1.5 10- 10 exp(-4640/T gas ) cm 3 s" 



(8) 



where the two terms are for the CHJ and CH + products, respec- 
tively. The C/C n ionisation balance then reads 

([C] - Xc)h = x c x e n H k r + x c xp A}i n H k x + x c yn H k a /2, (9) 

where y = 2n H Jn H is the H2 fraction and [C] = n c /n H = 
n(C i)/riH + xc is the total carbon abundance, which we take to 
be 1.3 X 10~ 4 . This leads to an expression for xc'. 



x c = [C] 1 + x, 



n H k r 



n H k x n H k a 



+ y- 



2k, 



(10) 



Combining Eqs.[3]and 10 yields a third-degree equation for x e 
whose solution is x s . In shielded regions, xc drops rapidly to 
form C and then CO. The electron fraction x e follows this evolu- 
tion down to a value Xj c corresponding to the case of dark clouds 
where the ionisation is mostly due to cosmic rays. In this case, 
we assume that the electron fraction follows the formula that 
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Fig. 1. Ratio of ion fractions calculated as decribed in Section 
[3] and of ion fractions calculated omitting the recombination of 
Cn with PAHs {k x = 0). Results are shown for n# = 0.1 (blue), 
30 (green), 100 (red), 1 000 (magenta), and 10000 ctrT 3 (black), 
for xh (dashed lines from bottom to top) and xq (solid lines from 
top to bottom). 



Williams et aL] ( |l998| l derived from an analysis of C ls O, H 13 CO + 
and DCO + observations in low-mass cores: 



x dc = 2000 



CR 



\ yn H 



(11) 



where y - 2n^ /tin- We then write x e - MAX(xd c , x s ). 

The Cn and Hn fractions are then derived from Eqs. [3] 
and [TO] Fig. [TJ shows the influence of the recombination of Cn 
with PAHs on the calculation of x H and x c . When it is omitted 
(k x = 0), the impact is the strongest when the radiation field in- 
tensity is the weakest. The C n fraction is over estimated by a 
factor of up to 20 in dense clouds, whereas the H 11 fraction is 
underestimated. The recombination with PAHs has no influence 
when Go > 50, whatever the density. 

The above fractions depend on the gas temperature, T gas , as 
well as on the hydrogen molecular fraction y. To derive these 
quantities, we performed a grid of simulations with CLOUDY 
(Ferl and et al.| 1998). We then interpolated on these grids to ob- 
tain r gas and y in the optically thin limit (Nh < 10 20 H/cm 2 ). We 
note that even if the gas is fully molecular (y = 1), the reactions 
of C 11 with H2 have little influence because their rates are much 
lower than those of reactions involving PAHs. The results are 
shown in Fig. §for a gas density of n H = 100 cm . 



4. Influence of local physical conditions 

4.1. Variations with the gas density 

We study the variations of the spinning dust spectrum for grains 
illuminated by the standard interstellar radiation field, ISRF, 
with Go = 1, and for clouds with densities ranging from n# = 
0.1 to 10 5 cm -3 . For each density, the gas properties are esti- 
mated according to n# and the radiation field intensity, as de- 
cribed in Section [3] The resulting spectra are shown in Fig. [3] 
and we also display the variation of the peak frequency and of 
the maximum intensity of the spectra as functions of n#. 

When tin < 10 cm 3 , the rotational excitation and damping 
are dominated by the IR emission (Fig. [4]). Thus the peak inten- 
sity of the spinning dust spectrum does not much depend on the 
gas properties and it varies little when increasing n#. The slight 




Fig. 2. Ion fractions xh (dashed line) and xq (dot-dashed line), 
H2 fraction y (dotted line), and gas temperature T gas (full line) as 
a function of the intensity of the radiation field Go for n# = 100 
H/cm- 3 . 
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Fig. 3. a) Spinning dust spectra per H column density for in- 
terstellar PAHs illuminated by the standard ISRF, Go = 1, for 
n H = 400, 1000,3 000,4000,7500,20000 and 50000 cm- 3 , b) 
Peak frequencies (full line + circles) and intensities per H col- 
umn density at the peak (dashed line + diamonds) of the spinning 
dust spectra for Go = 1 and n# = 0.1 - 10 5 cm- 3 . 



increase of the peak frequency and of the intensity of the spectra 
for tin > 1 cm 3 is caused by the increase of the influence of col- 
lisions with neutral species (Fig. If 10 < n# < 4000 cirr 3 , 
the gas-grain interactions become the dominant processes for 
rotational excitation and damping. This results in an increased 
peak frequency when h# increases. For 4000 < n# < 10000 
cm -3 , the gas-grain interactions are still the dominant processes 
but the H 11 and C 11 fractions dramatically decrease when n# in- 
creases. As a result, the ion-grain collisions are less efficient and 
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n H (cm-3) 

Fig. 4. Contribution of the different processes to the damping 
(solid line) and the excitation (dashed line) of the spinning dust 
emission for a grain with a radius of 3.5 A illuminated by the 
standard ISRF with Go = 1. The contributions are normalized to 
the drag produced by non-elastic collisions in a pure H i gas of 
a density n# (Draine & Lazarian[T9 98). The red lines show the 
contribution of IR emission, the green lines of collisions with 
neutrals, the blue lines of collisions with ions, the black lines 
of plasma drag, the magenta line of photoelectric effect, and the 
yellow line of H2 formation. 



the peak frequency decreases (Fig. [4]). Finally if > 10000 
cm" 3 , the peak frequency increases, while the emissivity varies 
only little. 

The numbers given here are valid only for the ISRF with 
Go = 1. If the intensity of the radiation field is increased (de- 
creased), the threshold values from one domain to the next are in- 
creased (decreased). We emphasize the importance of correctly 
estimating the ion fractions when modelling spinning dust emis- 
sion. 

4.2. Variations with the intensity of the interstellar radiation 
field 

We study the variations of the spinning dust spectrum for grains 
illuminated by the ISRF scaled by Go factors between 10~ 4 and 
10 4 . Several gas densities are considered, from 0.1 to 10 4 H/cm 3 . 
The ion fractions are calculated according to the values of h# and 
Go (see Section[3]l. The results can be seen in Fig. [3] 

First, for low Go the dominant processes are gas-grain in- 
teractions: the excitation is led by collisions with ions, while 
the damping is dominated by plasma drag and collisions with 
neutrals. This explains why the spinning dust spectrum varies 
very little with Go- Depending on the gas density, this applies to 
Go < 10~ 3 - 10. A correct estimate of ion fractions is thus im- 
portant. For instance, if the recombination of C 11 with PAHs is 
omitted when calculating Xh and xq (Fig. [TJ, the spinning dust 
spectrum is shifted to higher frequencies for low-Go- The peak 
shift caused by reaction 6 is between 1 and 25 GHz for n H — 30 
to 10000 cm" 3 , while the intensity is increased by 5 to 45%. 
A second domain can be distinguished for 1 < Go < 100. The 
grains become positively charged when Go increases. For these 
intermediate radiation field intensities the collisions with ions 
are still the dominant exciting process even if less efficient. As 
a result, we see a slight decrease of both the emissivity and the 
peak frequency of the spinning dust emission, which is more 
intense when the density is higher. For Go > 100, IR emission 




10 4 10 3 10 2 icr* io° 10 1 10 2 10 3 10 4 

Go 



Fig. 5. a) Peak frequencies of the spinning dust spectra for Go = 
10~ 4 to 10 4 and n H = 0.1 (blue), 30 (green), 100 (red), 1000 
(turquoise), and 10000 cirT 3 (purple), from bottom to top. b) 
Intensities per H column density of the spinning dust spectra at 
peak frequencies for the same environmental parameters. 



dominates the rotational excitation and damping, so that increas- 
ing Go results in increased emissivity and peak frequency. 

We stress that the variations of the spinning dust spectrum 
with Go are not monotonous and that they depend on the gas 
density. Moreover, the PAHs mid-IR emission is proportional to 
Go, so that we do not expect a tight correlation of the IR bands 
with the AME if the latter is caused by spinning dust emission. 

4.3. Influence of the cosmic-ray ionisation rate ( C r 

Gas-grain interactions through collisions or plasma drag are one 
of the dominant processes for the rotational damping and excita- 
tion. A proper evaluation of the ion fractions in the gas, xh and 
xc, is therefore important for the estimation of the spinning dust 
emission. The carbon is ionised by the UV radiation field and 
thus xc depends on Go (see Section[3]l. However, because of the 
Lyman cut of the radiation field, inside clouds the hydrogen is 
mainly ionised by cosmic -rays and x H depends on the cosmic- 
ray ionisation rate g. In the previous sections, we adopted the 
Galactic standard value, (cr - 5 x 10~ 17 s _1 , measured towards 
dense molecular clouds. However, several authors have inferred 
an enhanced cosmic -ray ionisation rate from observations of dif- 
fuse clouds {tin < 500 cm 3 ). They showed that this rate can 
exceed 50 times the standard cosmic-ray ionisation rate ( |McCall 
|et al.|2003[|Liszt|2003l|Siiaw et al.|2006> . 

To investigate the influence of the cosmic-ray ionisation rate, 
we modelled spinning dust spectra for n H — 30 cm 3 , Go = 1, 
and f = CRx £ C r with 1 < CR < 100 (Fig. |6). When ( CR 
is increased by two orders of magnitude, xh is multiplied by 
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n H = 30 cm" 
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Fig. 6. Peak frequencies (solid line + circles) and intensities per 
H column density at the peak frequencies (dashed line + dia- 
mond) of the spinning dust spectra as a function of the cosmic- 
ray ionisation rate equal to the standard rate £cr = 5 x 10~ 17 s _i 
multiplied by the factor CR = [1-100]. 

10 and the spinning dust spectrum is shifted by 10 GHz and its 
intensity is multiplied by three. The cosmic-ray ionisation rate 
therefore has a significant impact on the spinning dust emission. 
Its variations as a function of the gas density may produce a more 
complex behaviour than those presented in Figs.[3]and[5] 

5. Spinning dust emission with radiative transfer 

An interesting result obtained with the Planck data is that there 
is AME associated with all the interstellar phases of our Galaxy. 
It is detected in diffuse and dense media, ionised or neutral, 
and its spectrum can be fitted with a basic spinning dust model 
in all cases ( |Planck Co llaboration et al. 2 01 la|b) >. However, at 
finer angular resolution (~ a few arcminutes), the morphology 
of molecular clouds differs in the mid-IR and in the microwave 
( |Castellanos et al.pOTT) |Vidal et al.||20TTT >. We showed above 
that the spinning dust spectrum is very sensitive to the environ- 
mental conditions and we claim that it could explain the ob- 
served differences. Indeed, the variations of the gas properties 
and Go in dense molecular clouds are strong from the cloud edge 
to the centre. Hence, dense clouds are excellent targets to test the 
hypothesis of the spinning dust emission at the origin of the ob- 
served AME. Therefore, we investigate the spinning dust emis- 
sion emerging from dense interstellar clouds using the combina- 
tion of DustEM and CRT, in addition to SpDust. 

We model starless spherical clouds with density distributions 
that are suitable to fit interstellar clouds ( |Arzo umanian et al.| 
20TT]|Dapp & Basu|2009| l: 



n(r) 



no 



1 + (r/#o) 2 



for r < R n 



(12) 



with R oM the outer radius, Hq = R/3 the flat internal radius, and 
no the central density. The clouds are illuminated by the ISRF 
extinguished by a visual extinction Ay — 1. We consider that 
small and large grains are present everywhere in the clouds with 
constant abundances. The gas ionisation state is estimated as de- 
scribed in Section[3] 
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Fig. 7. Gas temperature as a function of the cloud radius for 
clouds with R ou , = 0.1 pc and no = 10 4 (solid line) and 10 5 
H/cm 3 (dashed line), and for a cloud with R oul — 1 pc and 
no = 10 3 H/cm 3 (dotted line). 

_S (MJy/sr) 
0.080 




Fig. 8. Surface brightness maps at 12 fim (a) and 30 GHz (b) and 
the emission profiles (c) at 12 jjm (solid line), 12 GHz (dotted 
line), 30 GHz (dashed line), and 44 GHz (dot-dashed line) for 
the starless molecular cloud with «o = 10 3 cirT 3 . 



to estimate r gas in spherical model clouds. This model takes into 
account cosmic -ray heating, line cooling, coupling between gas 
and dust, and photoelectric heating. The calculations follow the 
description of Goldsmith (2001 ) with the exception that the line 
cooling rates are estimated with Monte-Carlo radiative transfer 
modelling instead of using the large velocity gradient approx- 
imation ( jJuvela &~Y sard 2011 1. The line transfer is calculated 
with the abundances listed in Goldsmith (2001 ) and assuming a 



turbulent linewidth of 1 km/s. The resulting temperature profiles 
are shown in Fig. [7] 



5.1. Gas temperature 

In the dense molecular clouds the gas temperature grid calcu- 
lated with CLOUDY in the optically thin limit is no longer rel- 
evant (see Section [3]). We used our own radiative transfer code 



5.2. Starless molecular cloud 

We model a starless molecular cloud with an outer radius R out = 
1 pc and a central density no = 10 3 H/cm 3 . Fig. [8] shows the 
emerging dust emission for this cloud. Both the PAHs mid-IR 
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molecular cloud with no = 10 3 cm 3 (dotted line), 10 4 cm 
(dashed line), and 10 5 cm -3 (solid line). 



0.0 
R/Rout 



emission and the spinning dust emission increase towards the 
centre of the cloud. However, the shapes of the profiles are quite 
dissimilar: the mid-IR profile is approximately twice as wide as 
the microwave profile. In this cloud, the intensity of the radiation 
field varies from 0.24 at the edge to 0.06 at the centre. With these 
low intensities, the dominant processes for the excitation and the 
damping of the spinning motion of the grains are the interac- 
tions with the gas. The spinning dust emission is consequently 
stronger at the centre of the cloud. On the other hand, because 
the PAH mid-IR emission is directly proportional to the inten- 
sity of the radiation field and to the column density, the grains 
emit more at the edge of the cloud, which produces a broader 
emission profile. 



5.3. Dense molecular globules: spinning dust as a tracer of 
grain growth from diffuse to dense medium 

We modeled two dense molecular globules with an outer radius 
R = 0.1 pc and central densities no = 10 4 and 10 5 H/cm 3 . The 
emerging dust emission at 12 yum, at the centre of two Planck- 



LFI filters (30, 44 GHz), and at 12 GHz are shown in Fig. [10 
As expected, the brightness maps in the mid-IR and in the mi- 
crowave are dissimilar and even anti-correlated depending on the 
central density and on the frequency. The emission profiles in 
the microwave range have different widths depending on the fre- 
quency. This reflects the variation of the peak frequency of the 
spinning dust emission with the local physical conditions at a 
given radius (Fig.|9]l. Apart from these general trends, the results 
are strikingly different depending on the central density of the 
cloud. We emphasize that these findings are quite different from 
what is expected and observed at degree scale in the diffuse in- 
terstellar medium where the mid-IR emission and the AME are 
correlated. 

The cloud with «o = 10 4 H/cm 3 illustrates the complexity 
of the variations of the spinning dust emission with the environ- 
mental properties (Figs.[H)fc-c). The peak frequency of the spin- 
ning dust emission varies from ~ 45 GHz at the edge to ~ 33 
GHz at the centre of the cloud (Fig. [9j, while 0.1 < Go < 0.25. 
These peak frequencies ar e hig her than in the previous case with 
«o = 10 3 H/cm 3 (Section 5.2 1, whereas the mid-IR emission is 
partly extinguished towards the centre. This explains why the 
profiles at 30 and 44 GHz are wider than the mid-IR profile. The 




Fig. 10. Surface brightness maps at 12 //m (a) and 30 GHz (b) 
and the emission profiles (c) at 12 fim (solid line), 12 GHz (dot- 
ted line), 30 GHz (dashed line), and 44 GHz (dot-dashed line) 
for the dense molecular globule with «o 
lower figures (d, e, f) show the same for no = 10 5 cm" 3 



10 4 cm 3 . The three 



width of the 12 GHz profile varies only little compared to the 
previous case because of its distance from the peak frequency. 

For the densest cloud, with no = 10 5 H/cm 3 (Figs. lOi-f), 
the brightest region in the mid-IR is a ring at the edge of the 
cloud. The emission profile of these clouds can be explained in 
two ways. Either the PAHs are present throughout the cloud with 
a constant abundance but the radiation field is too extinguished 
at the centre for the grains to emit efficiently in the mid-IR. 
Some of the emitted IR photons also are absorbed by the dust, 
which additionally decreases the surface brightness at the centre. 
Alternatively the abundance of small grains decreases from the 
edge to the centre of the cloud. There is no way to distinguish 
these two scenarios with only data in the mid-IR. However, this 
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ambiguity is lifted with a microwave emission map. Indeed, if 
PAHs are present at the centre of the cloud, the rotational excita- 
tion caused by the interactions with the gas particles is sufficient 
for them to emit in the microwave range. The peak frequencies 
vary from ~ 30 GHz at the edge to ~ 27 GHz at the centre, which 
coincides with the Planck-LFI bands and explains why the emis- 
sion increases towards the cloud centre for the three frequen- 
cies. Because the dense clouds are transparent in the microwave 
range, spinning dust emission from the cloud centre can be ob- 
served. On the other hand, if small grains are absent from the 
centre, only a ring of mid-IR and spinning dust emission is ob- 
served at the edge of the cloud. As a result, spinning dust emis- 
sion could be used to trace the evolution of small grains from the 
diffuse to the dense medium because the comparison between 
mid-IR and microwave data allow us to trace the variations of 
the PAHs abundance. In this context, the potential of comparing 
the AME to the mid-IR bands has recently been illustrated by 
|Tibbs etal.| ( |20TT] i. 

6. Conclusion 

To understand the emerging spinning dust emission from inter- 
stellar clouds, we studied the influence of the gas ionisation state 
and of radiative transfer on the grain rotational motion with con- 
sistent models for the gas state and for the interstellar clouds. We 
found that the spinning dust spectrum is sensitive to the abun- 
dance of the major H II and C n ions and that these abundances 
must be described consistently in the modelling by including a 
treatment of the gas state. In addition, we showed that radia- 
tive transfer within interstellar clouds has surprising effects, e.g., 
that the mid-IR emission and the AME can be anticorrelated to- 
wards the cloud centre. From these findings, we argue that the 
AME from dense interstellar clouds provides new ways of grain 
growth from diffuse to dense medium (e.g., depletion of small 
grains by coagulation on larger ones), providing it is observed at 
relevant (a few arcminutes) angular scales using ground-based 
radio telescopes. 
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